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ABSTRACT 

Two atmospheric probing experiments 
involving beyond-the-horizon propagation of 
microwave signals are reported. In the 
first exneriment, Doppler-shif t caused by 
the cross-path wind is measured by a phase- 
lock receiver with the common volume dis- 
placed in azimuth from the great-circle. 
Variations in the measured Doppler shift 
values are explained in terms of variations 
in atmospheric structure. The second 
experiment makes use of the pseudorandom 
sounding signal used in a RAKE communication 
system. Both multipath delay and doppler 
shift are provided by the receiver permit- 
ting the cross section of the layer 
structure of the atmosphere to be deduced. 
The angular dependance and the cross-path 
wind in each layer are displayed. 


1. INTRODUCTION 

This paper is concerned with the problem of identifying the structure and 
dynamics of the atmosphere from the coherent detection and analysis of sound- 
ing signals after their transmission through the scatter medium. Section I 
discusses the coherent CW soundings performed on the UW-Collins 230 km link. 
This work, which relies on the antenna beams for spatial resolution, has 
successfully identified cross-path winds in scattering strata at different 
heights and has led to a model which hopefully improves the understanding of 
the role of cross-path winds and atmospheric layer structure in the propaga- 
tion mechanism. A model is presented which predicts the doppler shift of the 
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received signals in terms of the cross-path motion of the scatterers through 
the family of ellipsoids of constant path delay. This model also predicts a 
unique amplitude vs. path delay and doppler-shif t relationship for scatterers 
that are moving with a given horizontal cross-path wind speed in each strata. 

Successful tests of amplitude vs. doppler shift for this model were per- 
formed by synchronously pointing the UW and Collins antennas off the great 
circle and noting that the average doppler shift increased appropriately with 
antenna beam displacement from the great circle. 

A more useful test of the model which is discussed in Section II became 
possible with the results of certain RAKE channel sounding experiments describ- 
ed by Barrow, et al. , (1965) and Abraham, et al., (1967). Interpreting the RAKE 
data in terms of the model yielded layer heights and wind speeds that compared 
favorably with the weather bureau sonde data available at the nearest times and 
location. Since the RAKE technique achieves spacial resolution by the use of 
wide-band, short correlation-time sounding signals rather than by small antenna 
beams, several system advantages are realized. These are discussed in Section 
3, together with a hybrid system which would employ a vertical fanbeam. 

2. CW SOUNDING ON THE UW-COLLINS LINK 

In the UW experiment, a phase-stable 960 MHz, CW signal derived from a 
1 MHz standard is transmitted at 10 kw from a 28 ft. steerable paraboloid at 
the Collins Communication Research Facility (CCRF) near Cedar Rapids, Iowa, to 
the University of Wisconsin Facility (UW) near Arlington, Wisconsin (Figure 1) . 
The signal received at UW on a 28 ft. steerable paraboloid is phase-coherently 
retransmitted at 810 MHz and 1 kw from the same (duplexed) antenna to CCRF. 

The transponded signal is received at CCRF on a separate 28 ft. paraboloid 
slaved to the CCRF transmitting antenna. The round-trip signal is phase com- 
pared with the signal from the original standard in a linear phase-detection 
system. The phase difference, due almost entirely to the propagation path, is 
continuously recorded on two simultaneously available ranges ( 17 . 12 tt and 60 tt 
radians). Each phase output automatically resets at its range limits, thereby 
providing linear, piece-wise continuous data over an unlimited number of 
cycles, (Figure 2). Residual phase noise of the closed-loop system is about 
30 rms. Experimental error introduced by the antenna pointing systems is con- 
sidered negligible. Antenna pattern measurements revealed a constant great- 
circle misalignment of approximately 0.2°. Data illustrated in Figures 3 
through 9 have not been adjusted for this small error. 

The signal phase behavior includes a wide range of variations 
(Birkemeier, ejt. al. , 1965). Certain features are consistently noted. Phase 
jumps 1, of up to, but not exceeding, tt radians occur coincident with deep 
fades, indicating the interference nature of the fading phenomena. With the 
exception of these phase jumps and occasional periods of aircraft interference, 
the instantaneous phase rate never exceeds that which can be produced by tropo- 
spheric winds. When the antennas are aligned on the great-circle azimuth the 
phase has essentially zero average slope. Random variations over several 
minutes rarely exceed about 10 tt radians (Figure 2) . When the antennas are 
synchronously set to positions off the great-circle azimuth, the phase accumu- 
lates continually in one direction (Figure 2) . Random phase variations about 
the average slope also become more rapid, suggesting that the phase spectrum 
of the received signal is dependent on the antenna pointing angle. 

The average slope of the phase-versus-time record corresponds to an aver- 
age doppler shift and at a given antenna angle is clearly defined over a few 
minutes (Figure 2) . 

A series of beam-swinging experiments revealed the average doppler shift 
to be a systematic function of the antenna pointing angle and the direction 
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Table 


I 


Beam Swinging Experiments 

Time at 


Exper intent 


Oate 


Time (CST) 

Each Anqle* 

A 

10 

June 

1966 

1437-1750 

1 

B 

7 

July 

1966 

1015-1430 

2 

C 

21 

Sept . 

1966 

1211-1313 

1 

D 

28 

Sept . 

1966 

0431-0626 

3 

E 

28 

Sept . 

1966 

1840-1915 

4 

F 

29 

Sept . 

1966 • 

0230-0430 

3 

G 

29 

Sept . 

1966 

1055-1145 

5 

H 

24 

July 

1967 

2225-0230 

5 

1 

25 

July 

1967 

0440-0710 

6 

J 

25 

July 

1967 

1010-1 1 15 

7 

K 

25 

July 

1967 

1115-1200 

7 

L 

25 

July 

1967 

1200-1315 

7 


Sequence 
of Angles— 

1 

2 

3 

3 

3 

3 

4 
3 

5 
3 
3 
3 


*Time at Each Angle 

1 . 10 minutes . 

2. 15 minutes. 

3. 10 minutes between 0.6° N & S, 5 minutes each larger angle. 

4 . 1.5 minutes . 

5. 10 minutes between 0.4° N & S , 5 minutes each larger angle. 

6. 4 minutes. 

7- 3 minutes. 


^Sequence of Angles 

1. 0.2° HIS, 0.4° N 5. S. etc rt 

2. 0.4° HIS, .1.0° N ( S, 1.0 NtS, 1.0° N 1 S. 0.6° NSS. 

3 . Sequential ly S to N. 

4. Sequentially N to S. 

5. Sequentially S to N from 5.0° S, then 5.0° S to 7 .0° $. 


TIME (CST) 

1500 1600 1700 1800 1900 



Figure 4. Measured average Doppler 
shifts vs. azimuthal antenna-pointing 
angle and time. 

and magnitude of the crosspath wind. In each experiment the antennas were syn- 
chronously set for several minutes at a sequence of azimuthal angles. Experi- 
mental observations are summarized in Table 1 and Figures 5 through 9. 

Doppler shifts measured at equal angles on opposite sides of the great 
circle are essentially Identical in magnitude in Figures 5, 6 and 8. The sign 
of observed doppler shifts agrees with the direction of crosspath winds in the 
lower portion of the common volume, i.e., the received frequency is greater 
(less) than the transmitted when the antennas are pointed upwind (downwind) . 

The symmetric, systematic, non-linear characteristic of these doppler shift 
versus antenna pointing angle curves cannot be ascribed to nonstationary mete- 
orological processes . The daily to hourly variation in this curvature suggests 
dependence on the prevailing reflective angular dependence of the atmospheric 
refractive structure. 

The observed data will be discussed in the context of the proposed doppler 
shift model which is developed in terms of the average crosspath wind. Concur- 
rent wind data were obtained by tracking meteorological pilot balloons with a 
theodolite at launch intervals of 30 minutes to one hour. The analysis of 
tracking data included computation of vertical profiles of crosspath wind 
speed. Representative crosspath vertical profiles are shown in Figure 10. 
Balloons were normally launched from a field site at midpath beneath the com- 
mon volume. 

2.1 Doppler Shift Model 

To each point in space there corresponds a ray-path length L(x,y,z) which 
is defined by a ray from the transmitter to receiver passing through the given 
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point. Surfaces of constant path length or phase are approximately ellipsoids 
of revolution having the radio terminals as foci (Figure 11). A signal obtain- 
ed via any single moving point is shifted in frequency by an amount proportion- 
al to the rate at which the point cuts equiphase surfaces. Thus the "single- 
scatterer" doppler shift f is given by 

f = - -f %: = - f V • VL (1) 

A dt A 

where X is the wavelength and V is the velocity of the moving scatterer. As- 
suming the simplest ray path of two straight lines through the scatterer, L is 
determined using the coordinate system illustrated (Figure 11) with origin at 
the midpoint of the straight (chordal) line between the terminals. 


L = [ (d + x) 2 + y 2 

+ z 2 ] 1 / 2 + t(d - x) 2 + y 2 + z 2 ] 1 / 2 

(2) 

Applying (1) to (2) yields 



f - A_ 

ux(l - 4d 2 /L 2 ) + vy + wz 

(3) 

f XL 

1 - 16x z d z /L H 


where (u,v,w) are the (x,y,z) components of V. 

The relative importance of the various terms in (3) depends on the posi- 
tion of the scatterer and on the relative magnitude of the velocity components. 
The along-path (u) contribution is negligible because the ellipsoids are longi- 
tudinally very "flat". The factor (1 - 4d 2 /L 2 ) is small, and equal to 4z|,/L 2 = 
4h^/d 2 , where z m and h ra are the midpath heights of an ellipsoid above the 
chordal axis and spherical earth, respectively. This factor ranges from 10 -t * 
to 10 -2 from the bottom to the top of the common volume for the path employed 
in these experiments. 

At a given point, the average horizontal velocity normally greatly 
exceeds the average vertical velocity. We interpret our observations in terms 
of drift of scatterers with the average crosspath velocity v. Under these con- 
ditions, (3) predicts for a single scatterer in the vicinity of midpath, 

f = - — v (— £■) = - — v sxn a = - — va (4) 

where a is the azimuthal angle of the single scatterer from the great circle. 
For a given crosspath wind speed , pointing the antennas to azimuthal angles 
away from the great circle tends to favor scatterers at larger angles a which 
consequently increases the observed doppler shifts. There are many scatterers 
at various positions, however, which create a spectrum of doppler shifts in the 
received signal. It is necessary to specify the doppler spectrum and its de- 
pendence upon azimuthal antenna pointing angle a a in order to explain the shape 
of the curves of average doppler shift f vs. a a . 

Consider numerous scatterers moving with the wind, each contributing a 
signal component having a doppler shift given by (4). Because of the longitu- 
dinal flatness of the ellipsoids, scatterers even widely displaced from midpath 
may be included with little error. For a given constant crosspath speed the 
doppler shift varies linearly with crosspath position and is independent of 
height. The shape of the doppler spectrum which results is principally deter- 
mined by the relative strengths of scattered signals as a function of crosspath 
position. The signal power from a given position depends upon the combined an- 
tenna illumination pattern and upon the scattering cross section at the posi- 
tion. Let G(a) denote the multiplied azimuthal pattern function of the anten- 
nas and W(a) the scattered power per unit azimuthal angle. For synchronously 
pointed antennas with boresights at a a . 
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Figure 5. Measured average Doppler 
shifts vs. azimuthal antenna-pointing 
angle. 



Figure 7 . Measured average Doppler 
shifts vs. azimuthal antenna-pointing 
angle. 



Figure 6. Measured average Doppler 
shifts vs. azimuthal antenna-pointing 
angle. 



Figure 8. Measured average Doppler 
shifts vs. azimuthal antenna-pointing 
angle. 
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S (f ) df - S a (a,a a ) da = G(a-a a ) W(a) da (5) 

where S(f) is the doppler spectrum S a (a,a a ) is the doppler spectrum written in 
terms of a, and f and a are related by the scaling relation (4). 

Bello (1965) gives the average doppler shift f in terms of the doppler 
spectrum as 

f = / f S(f) df/ / S(f) df (6) 

which according to (4) and (5) may be written 

f = - ~ f aG(a-a a ) W(a) da/ / G(a-a a ) W(a) da (7) 

This equation provides a relation between f and a a in terms of the antenna 
patterns and the "angular dependence" function W. W includes the azimuthal de- 
pendence of scattering cross section and/or reflection coefficient depending on 
the degree of anisotropy of the atmospheric structure involved. W is assumed 
to decrease monotonically and symmetrically with increasing angle from great 
circle. The doppler spectrum S is weighted by W toward the great circle from 
a a , and the average doppler shift is reduced from the value associated with 
scatterers at a a . 

An angular dependence function for reflectivity based upon the Tatarski 
(1961) model for a refractively turbulent medium is 

n = const sin 2 B / sin 11 ^ 3 (0/2) (11) 

where 0 is the scattering angle and 6 is the angle between the electric field 
vector and the ray to the receiver. The combined wind and phase data suggest- 
ed that dominant signal contributions were frequently propagated via a low lay- 
er of limited vertical extent. Supporting observations of tropospheric layer 
structure by Atlas, et . al . (1966) and Hardy, et . al . (1966) using radar back- 
scatter, and by Lane (1966) using radar and radio refractometers , stimulated 
analysis of a thin, turbulent scattering layer. Using n (with sin 2 f? = 1) and 
gaussian antenna patterns, the "doppler spectra" of Figure 12 were computed 
from (5). The peaks of S a (a,a a ) occur to the great-circle side of the pointing 
angle. These spectra are converted to doppler spectra by means of the scaling 
relation (4) for a given average crosspath wind velocity. Numerical integra- 
tion of (7) for such doppler spectra yields the set of average doppler shift 
versus antenna pointing angle curves of Figure 13 for several heights of the 
thin, turbulent scattering layer. These curves are symmetric in a a , a property 
of (7) for symmetric W and G. A change in the magnitude or direction of the 
crosspath wind changes the vertical scale or the sign of the curves but not the 
shapes. The foregoing doppler shift model provides the qualitative features 
and correct orders of magnitude required to interpret the data. 

2 . 2 Interpretation 

The crosspath wind profiles have in several cases permitted identification 
of the height range of significant scattering. 

During experiments conducted on 7 July 1966, 21 September 1966, 29 Septem- 
ber 1966 and 24, 25 July 1967, the crosspath component of the wind was from the 
northwest (NW) at all altitudes. The sign of the average doppler shift corre- 
sponded to this direction, i.e., was positive (frequency increase) for antennas 
directed toward the windward side of the great circle, and negative (frequency 
decrease) for the leeward side. During the experiments of 25 July 1967 in 
which the beams were azimuthally swung in both an unelevated and elevated posi- 
tion, the average doppler shift was greater for the cases of elevated beams 
(Figure 9). The observed increase was in agreement with greater wind speeds 
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ANTENNA POINTING ANGLE (DEG.) 

Figure 9 . Measured average Doppler 
shifts vs. azimuthal antenna-pointing 
angle. 



Figure 10. Typical crosspath wind pro- 
files measured during varied meteoro- 
logical conditions. 


Z(Km) 



Figure 11. Longitudinal and transverse 
sections of propagation link and of 
ellipsoids of constant path length. 

Path length Increment between ellip- 
soids is 80 X shown for X = 0.37 m. 



Figure 12. Theoretical Doppler spectra 
(written in terms of a) for different 
a a for a thin, turbulent scattering 
layer at 785 m above the earth, comput- 
ed using Tatarski scattering cross sec- 
tion and Gaussian antenna patterns 
(each with 3.0° half-power beamwidth) . 
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ANTENNA POINTING ANGLE a 0 , (DEG.) 


Figure 13. Theoretical average Doppler 
shift vs. antenna-pointing angle for the 
two-way path for thin, turbulent scat- 
tering layers, computed using Doppler 
spectra such as those of Fig. 12. Indi- 
cated heights are approximately at the 
bottom, middle and top of the nominal 
common volume, v = 1 m/s. 


Figure 14. Comparison of data of ex- 
periments K and L (Fig. 9) with theo- 
retical average Doppler shift for sev- 
eral layer heights. Wind speeds are 
average of those measured at the indi- 
cated heights during the experiment. 
Data are compensated for 0.2° constant 
misalignment . 


measured at the higher altitudes. 

Experiments on 28 September 1966 were conducted when the crosspath compo- 
nent was from the NW at all heights above approximately 400 m altitude but from 
the southeast (SE) between the surface and the reversal height (Figure 10) . 
During the experiments the sign of the doppler shift corresponded to NW winds. 
Thus although a large percentage of the total propagation path is at altitudes 
below 400 m, this region did not appear to contribute the significant doppler- 
shifted signal components. 

During the experiment of 10 June 1966, the crosspath component was SE be- 
low, and NW at all heights above 1.6 km. The sign of the average doppler shift 
corresponded to the SE winds in the lower portions of the common volume. 

A vertical beam-swinging experiment in which the crosspath wind reversed 
direction from SW to NW at a height of approximately 1.5 km was conducted on 
12 May 1967 and phase data is shown in Figure 15. With the antennas se.t 2.0° 
off the great-circle azimith and at 0.0° elevation, the sign of the average 
doppler shift corresponded to the SE winds in the lower zone. When the beams 
were elevated to 2.0°, the average doppler shift reversed to the sense corre- 
sponding to the NW winds above 1.5 km. This result suggests that a crossed 
pair of bistatic, doppler-sensitive systems, measuring both crosspath wind com- 
ponents, might be used to remotely measure winds. 
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In their entirety these observations demonstrate that the dominant signal 
contributions were reradiated in the lower regions of the common volume, as 
distinct from either the upper portions of the common volume or the regions be- 
tween the antennas and the common volume. This interpretation agrees with the 
conclusions of earlier studies based on statistical analysis of amplitude data 
by Ikegami (1960) and Fengler (1964). The next section on the RAKE data also 
confirms that stronger reflecting layers tend to be the lowest in elevation. 

In Figure 14, theoretical curves, similar to those of Figure 13 but com- 
puted for single layers at different heights, are compared with the data of 
Figure 9. For clarity, only two data curves are included. The theoretical 
curves are scaled for the measured winds at the labeled heights. At low alti- 
tudes, where the wind speed increased with height, the height region of agree- 
ment is well resolved. The height of agreement for the case of elevated beams 
is poorly defined because the wind velocity was essentially constant between 2 
and 3 km altitude. Good resolution was also obtained when the same procedure 
was applied to the data of Figure 6. The heights of agreement were about 1.5 
and 2.0 km for experiments D and E respectively. 

The model based on the Tatarski hypothesis is inadequate to explain fea- 
tures of the data of Figure 8 v?hich were obtained during a stable, stratified 
nocturnal atmosphere. Stable layers may produce highly anisotropic space cor- 
relation functions of index of refraction that yield spectral energy only at 
wave numbers nearly normal to the earth's surface. This results in quasi- 
specular components arriving only via paths near the great circle plane. For 
this reason these components have nearly zero doppler shift independent of 
their scatterer's crosspath speed. Their contribution to the doppler spectrum 
is accounted for by emphasizing the angular dependence function W at small a. 
Since many strata of scatterers with differing angular dependence functions and 
differing wind speeds would normally be included within the UW link's antenna 
beams, it is hazardous to try to associate the W function in a given experiment 
with a single layer. However, in the next section it will be shown that such 
an identification may be possible with RAKE probing even with large beams. 

3. RAKE TROPOSCATTER SOUNDING 

The experiments of Barrow (1965) and Abraham (1967) analyzed in this sec- 
tion were the first reported tests of the RAKE adaptive communication system 
applied to the troposcatter channel. 

In its normal communicating mode of operation this system employs an 
estimator-correlator receiver containing a stored alphabet of reference sig- 
nals. The decision as to which signal from the alphabet was transmitted in 
each signaling interval is determined by cross-correlating the received signal 
with each stored reference signal and selecting the one with the greater cor- 
relator output. 

Since the received signal is distorted by the channel, the correlation 
process can be enhanced by predistorting the reference signals in the alphabet 
in an identical manner. To allow this, the receiver must learn the channel 
transfer characteristic and implement it as a delay-line model into which the 
reference signals are fed before cross-correlation with the arriving signal. 

It is the data from the channel sounding operation that is interpreted here in 
terms of the atmospheric layer structure and the crosspath winds. 

To enable the RAKE receiver to learn the channel's transfer characteris- 
tic, i.e., its various paths, their delays, gains and doppler shifts, a 10 MHz 
wide phase-stable, periodic, pseudo-random sounding signal is transmitted con- 
tinuously. Cross-correlation of this received signal with its replica at the 
receiver provides a correlation function which is proportional to the received 
signal intensity at those correlator delays which match the various path delays 
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of the channel. This cross-correlation function is the path-delay profile. A 
delay resolution of 0.1 usee, is provided in these data. Furthermore, the dop- 
pler shift associated with each path shows up as the beat-frequency of the 
cross-correlator output at that path delay setting. Figure 16 from Barrow 
(1965) shows an example of such a cross-correlation function produced by lin- 
early increasing the reference signal delay at the rate of .1 usee, per second. 
Note the total channel path-delay profile width is about 2.8 yseconds. 

A systematic increase in doppler frequency with delay is clearly visible 
in the profile. This was plausibly interpreted by Barrow in terms of vertical 
wind velocities through the common volume, with higher vertical wind speeds at 
higher altitudes. However, it can be shown that the profile agrees closely 
with the result expected from a thin layer of scatterers moving essentially 
only with the horizontal wind. To show this, the doppler cycles in the profile 
were counted, starting at the first arrival, and the resulting function compar- 
ed to that predicted by the model. To get the predicted function, consider a 
scattering point at (y,z). The path delay at the point relative to (0,z) is 
given by (2) and amounts to 

t=(L-L.)/c= y 2 / c(d 2 + z 2 ) 1/2 = ay 2 (12) 

min 

Next, let the reference signal delay be varied linearly in real time t at 
the rate of p sec/sec. If the total number of doppler cycles appearing in the 
swept multipath delay profile is counted as a function of sweep time T, begin- 
ning with the first-arriving signals, the result is given by 

N(T) = f d (t) dt (13) 

where f (t) is the doppler beat frequency as a function of time. 

Recalling that for the model the doppler shift varies linearly with cross- 
path displacement y and that the delay time t varies as the square of y , we 
have from (4) and the relation t = pt , 

f (t) = bvy = bv(t/a)^ 2 = bv(p/a)^ t ^ (14) 

d 


where b, v, a and p are constants. 
Finally , T 

N(T) = K f Q 



dt = const 



(15) 


Figure 17 shows a log— log plot of N vs. T using the raw, counted data of the 
profile. A clear power-law dependence is demonstrated, with the empirical pow- 
er equal to 1.58. 

Further consideration of this multipath profile revealed the likelihood 
that the antennas of the RAKE link were both displaced about one— half beam- 
width to one side of the great-circle path. This was evidenced by the fact 
that, if the model is correct, the envelope of the multipath profile appears to 
violate the possible laws of angular dependence. The envelope peaks at about 
.8 pseconds and then falls with increasing delay. No AGC was used. The total 
envelope amplitude change is also about 6 db , similar to the combined antenna 
patterns, suggesting the profile is due to a horizontal scattering layer lying 
near the bore-site intersection of the beams. Using the known antenna eleva- 
tion angle, the suspected antenna position is shown superimposed on the cross- 
section of the ellipsoids of constant time-delay in Figure 18. 

Meanwhile, continued data analysis by the original authors developed the 
doppler spectra for the signals arriving at 0.1 sec. time-delay intervals. 
These spectra revealed a consistent negative doppler shift. See, for example, 
Figure 19a. The prevailing cross-path winds were determined from weather bu- 
reau data to be from the West, indicating azimuthal antenna displacement to 
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Figure 15. 
-23 tt rad /min 




Figure 16. quadrature conponents of swept nultipath profile from Barrow, 
1965. This profile followed scattering function 151 of Figure 19 
in tiire by 6 hours. 


COmpon( f nts of swe P fc multipath profile from Barrow, 1965, 
This profile followed scattering function 151 of Figure 19 in time by 6 hours. 
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RELATIVE TOTAL PROPAGATION DELAY (fiSEC) 

Figure 17. Plot of accumulated doppler-beat cycles vs. sweep time. 

Koclital Aatcnna Pattern 



Figure 18. Apparent antenna position for RAKE link. Also shown is the con- 
struction of a doppler-spectrum from a single O.ly sec time shell and a single 
scattering layer. 


349 


IDENTIFICATION OF STRUCTURE BY MICROWAVE SOUNDING 


the East of the path. This predicted antenna displacement was consequently in- 
vestigated by the original authors and the displacement to the East has been 
confirmed. The exact position is not known to this writer, however. 

3.1 Atmospheric Structure From RAKE Profiles 

In each of the RAKE scattering diagrams (Figures 19a, 23, 25) ridges of 
peaked doppler-spectra appear to follow parabolic loci in the time-frequency 
plane. Such loci are predicted by the model for each scattering strata by the 
fact that horizontally drifting scatterers increase their doppler shift linear- 
ly with y while their signal delay increases with y 2 . Combining (12) with (4) 
and using the dimensions of the 480 km RAKE path, we obtain the delay-doppler 
locus for point scatterers at a height z and moving with cross-path wind speed 
v as 

x = 21.8 f 2 / v 2 + z 2 /72 - .64 (16) 

Here t is in microseconds relative to the ellipsoidal shell at the grazing ray 
intersection which represents the scattering location of the first possible ar- 
riving signals; f is in Hz, v is in meters/second, and z is in kilometers above 
the link chord. 

To facilitate the interpretation of the RAKE profiles in terms of the im- 
plied layer structure, a program was written on the basis of the model to de- 
termine the theoretical ridge for an arbitrary layer height, thickness, and 
cross-path wind speed and angular dependence function. In this program the 
relative doppler power spectral density for a given delay shell is proportional 
to the vertical thickness of the scattering layer in the shell at a particular 
value of y (or f) , multiplied by the displaced antenna function and the angular 
dependence. Figure 18 shows the sketch of the spectrum obtained from a single 
shell. A theoretical ridge is shown in Figure 20 for a layer with Tatarski an- 
gular dependence, a height of 12.8 km, a thickness of 0.4 km, and a cross-path 
windspeed of 45 m/second. These values produce good agreement to the dominent 
ridge in the RAKE scattering function in Figure 19a. The spectrum at the ori- 
gin in the RAKE diagram is assumed to have originated from signals scattered 
from near the great-circle plane at the grazing ray height. This determines 
the reference delay for each higher layer. 

A number of additional parabolic tracks in this profile have been identi- 
fied by plotting local spectra maxima for each shell. (See Figure 19b). Fur- 
thermore, several sharp spikes along the zero-doppler shift axis imply the ex- 
istence of several specular layers as well. In an effort to compare these 
"radio identified" layers with the atmospheric structure in the vicinity of the 
link we resorted to the nearest weather bureau radio sonde profiles of wind, 
temperature, and humidity for the dates of the RAKE tests. Analysis of these 
data led to a set of probable layer heights together with their respective 
cross-path wind speeds. Using (16) each meteorologically identified layer's 
parabolic track was plotted and compared with the radio-deduced layer heights 
and winds. Figure 21 shows surprising agreement with the radio-detected lay- 
ers since the sonde data precedes the radio run by eight hours. The result is 
nevertheless intriguing. Figure 22 based on the same profile shows the radio 
layers drawn in cross section across the antenna pattern. Wind profiles from 
the three weather stations are included for comparison to the radio-derived 
wind in each layer. The dashed lines indicate radio layers for which only the 
zero-doppler shift point on its track could be identified. The fact that such 
layers are likely to be specular agrees with the stable state of the atmosphere 
at their location as indicated by the sonde data. 

In Figure 24 the radio layers corresponding to the parabolic tracks of the 
scattering function of Figure 23 are shown. It is apparent that the first 0.8 
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Figure 19b. Parabolic tracks of scattering diagram No. 151. 
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Figure 20. 


Figure 21. 



Theoretical scattering diagram for a layer of thickness = 0.4 km, 
height above chord = 12. 5 km, and cross-path wind speed = 45 m/s. 

-UOPPLLR FRLnUENCY (llz) 



Comparison of RAKE parabolic tracks with tracks derived from 
meterological sounding. Washington, D. C. sonde preceeded 
radio run No. 151 by 8 hours. 
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y (km from grcat-circlc) 

Figure 24. Atmospheric cross-section derived from Figure 23 (Scattering 
diagram 134). Radio data 2/17/65 2226 GMT, Sonde data 
2/18/65 



Figure 25. Scattering diagram No. 181 (Abraham, 1967). 
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microseconds are not included in the display. As such, the solid line segments 
indicate that part of the layer which contributed to the segment of the ridge 
in the scattering function. 

The scattering function of Figure 25 appears to have coincided with a 
frontal passage. The cross-section derived is shown in Figure 26. 

3.2 Considerations of Angular Dependence 

Most of the parabolic ridges analyzed so far show general agreement with 
the -11/3 power law. It is difficult to determine the power law exactly, how- 
ever, since neither the antenna pattern nor its position is known accurately. 
Further RAKE studies performed carefully, however, may well allow the angular 
dependence to be determined for particular turbulent layers. 

In certain RAKE data, strong angular dependence is implied by "spikes" 
along the zero-doppler shift axis in the RAKE scattering functions. An example 
is the scattering function shown in Figure 27. The analysis of this and simi- 
lar profiles is in progress to evaluate the exact angular dependence function 
and, if possible, the degree of anisotropy of the space correlation function of 
the index of refraction in the layer. 

3.3 Tilted Layers 

Under certain atmospheric conditions, scattering layers may be uniformly 
slanted with respect to the horizon or layers may contain uniformly slanted 
"scatterers" . If such scattering surfaces are visualized against the back- 
ground of the ellipsoids of constant phase of path delay, it is apparent that 
maximum reflectivity will tend to occur off the great-circle at the point of 
greatest tangency of the scattering surfaces to the ellipsoids. A RAKE para- 
bolic profile for such a layer (for antennas aimed along the great circle) 
should then peak to the side of the zero-doppler axis, indicating "slanted an- 
isotropy" in the space correlation function of the refractive structure. There 
is the possibility that the antennas of the RAKE link were aligned under such a 
slanted atmospheric structure, accounting for their fixed misalignment. 

4. PROPOSED HYBRID SYSTEMS 

4.1 Measurement of Vertical Wind Velocities 

An important limitation of the simplified model presented here, in spite 
of the apparent agreement with this particular RAKE data, is that it ignores 
turbulent velocity fluctuations in the wind. Such random variations in veloc- 
ity become quite important if the cross-path wind component is small compared 
with the total wind speed. The effect of the velocity fluctuations is to 
spread the doppler spectra at the various RAKE delays. This greatly compli- 
cates the problem of calculating layer thickness from the RAKE data. A Scheme 
to measure the vertical velocity fluctuations by eliminating doppler shift due 
to horizontal motion is suggested by the link geometry. Both spatial resolu- 
tion from the beams and time resolution from the RAKE are required. If the an- 
tennas are designed to produce vertical fan-beams sufficiently narrow to illu- 
minate only the first Fresnel zone along the great-circle plane, then according 
to (3), doppler shift can occur only from vertical scatterer motion. (Tests 
with beams of 0.3° width will be attempted on the UW-Collins link in the near 
future) Combining such fan-beams with a RAKE system to provide height resolu- 
tion should allow doppler-sensing of the vertical wind velocity field. Hori- 
zontal turbulent velocities may be measured similarly by azimuthally displacing 
the beams off the great circle, correcting for the vertical velocity fluctua- 
tions in the measured doppler shift fluctuations. 

4.2 Total Wind Measurement by Crossed Links 
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Figure 26. Atmospheric cross-section derived from Figure 25.. (Scattering 
diagram 181). Radio data 2/19/65, 1831 GMT, Sonde data 
2/19/65, 1200 GMT. 
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Figure 27. Scattering diagram No. 142 showing layer with strong angular 
dependence. 
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If two scatter links are arranged in parallel and spaced so that their el- 
lipsoids of constant phase intersect orthogonally at approximately the height 
of the grazing ray intersection, it appears possible to measure the vertical 
and horizontal wind velocities simultaneously at a common point midway between 
the paths. This measurement is accomplished simply by forming the sum and dif- 
ference of the phase outputs of the two links operating coherently on carrier- 
wave signals. Vertical profiles would be achieved by vertical beam swinging. 

5. CONCLUSIONS 

Forward-scatter atmospheric probing has been demonstrated to be potential- 
ly capable of providing useful structural information on the atmosphere. Much 
work needs to be done, but the combination of forward scatter’s great sensitiv- 
ity and the possibility of achieving spatial resolution through the use of sig- 
nal processing techniques rather than by the use of large antennas makes the 
possibilities exciting indeed. 
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